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Abstract. We review the scaling relations for the critical current density (Jc) in 
NbsSn wires and include recent findings on the variation of the upper critical field {Hc2) 
with temperature (T) and A15 composition. Measurements of Hc2{T) in inevitably 
inhomogeneous wires, as well as analysis of literature results, have shown that all 
available Hc2{T) data can be accurately described by a single relation from the 
microscopic theory. This relation also holds for inhomogeneity averaged, effective, 
i7*2(r) results and can be approximated by i?c2(i)/^c2(0) = 1 - with t = T/T^. 
Knowing H*2{T) implies that also Jc{T) is known. We highlight deficiencies in 
the Summers/Ekin relations, which are not able to account for the correct Jc{T) 
dependence. Available Jc{H) results indicate that the magnetic field dependence for 
all wires from fioH 1 T up to about 80% of the maximum Hc2 can be described 
with Kramer's flux shear model, if non-linearities in Kramer plots when approaching 
the maximum Hc2 are attributed to A15 inhomogeneities. The strain (e) dependence 
is introduced through a temperature and strain dependent H*2 (T, e) and Ginzburg- 
Landau parameter ki(T, e) and a strain dependent critical temperature T'c(e). This is 
more consistent than the usual Ekin unification of strain and temperature dependence, 
which uses two separate and different dependencies on H*2{T) and H*2{e). Using a 
correct temperature dependence and accounting for the A15 inhomogeneities leads to 
the remarkable simple relation J^{H,T,t) = {C / ^ioH)s{e){l-t^-^^){l-t^)h^-^{l-hf , 
where C is a constant, s(e) represents the normalized strain dependence of i?*2(0) and 
h — H/H*2{T, e). Finally, a new relation for s(e) is proposed, which is an asymmetric 
version of our earlier deviatoric strain model and based on the first, second and third 
strain invariants. The new scaling relation solves a number of much debated issues with 
respect to Jc scaling in NbsSn and is therefore of importance to the applied community, 
who use scaling relations to analyze magnet performance from wire results. 
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1. Introduction 

Critical current measurements of NbaSn wires are commonly interpolated and 
extrapolated using a set of empirical relations for Jc{H,T,e), generally referred to 
as Summers scaling [1]. Although these relations show acceptable accuracy for the 
description of limited data ranges [2], large discrepancies appear between the fitted 
values for Tc{0) and i?c2(0) and actual measurements [3]. Also at temperatures 
above about 13 K and for large compressive and tensile axial strains the Summers 
scaling is inaccurate. It is thus not suitable for precise extrapolations of measured 
results across a wide range. These deficiencies are generally recognized and lead to 
continuing discussions in the literature on improvements [2-12]. Presently, much of 
the interpretation of measured results and ensuing discussions remain largely empirical 
and the proposed relations depend, for a significant part, on the initial assumptions. 
A notable exception is the recent work by Oh and Kim [11, 12] who use Eliashberg 
theory [13, 14] in their attempt to arrive at a better founded scaling formalism. 

As an alternative to empirical approaches, a series of relations resulting from 
microscopic theory exist for Type II superconductivity. Although the critical current 
cannot be related quantitatively to the microstructurc, various bulk pinning force 
models have been developed [15-23]. These can be combined with theory on strain 
dependence [24-26] and the field-temperature phase boundary [27-47]. The latter have 
been shown to reasonably describe measurements on well defined, quasi-homogeneous 
laboratory samples [35,43-46,48-51]. 

It is questionable, however, whether such relations can directly describe the 
behavior of NbaSn in technical wires due to the inevitable presence of inhomogeneities. 
Compositional inhomogeneities result from the A15 stability range from about 
18 to 25 at.% Sn [52], combined with the solid state diffusion processes during 
the A15 formation reaction in wires. The resulting Sn gradients have been 
unambiguously detected in wires by compositional analysis [53-57], property gradient 
measurements [54-56, 58-60] and simulations [56, 60-62]. These compositional 
inhomogeneities lead to significant changes in the electron-phonon interaction strength 
Aep [63]. The interaction strength changes from the weak coupling BCS limit (Aep <S 1) 
for 18 at.% Sn to Aep = 1.8 for the stoichiometric composition [48,63]. As a result, 
significant changes occur in Hc2{T) with Sn concentration [48,49,63-67]. Since all 
stable A15 compositions from 18 to 25 at.% Sn are observed in wires, a distribution 
of properties of about 5 T < fiQHc2{Q) < 31 T and 6 K < Tc(0) < 18 K will occur 
over the A15 volumes [56]. In addition, it can be expected that also the strain state 
is far from homogeneous over the A15 volumes in a wire. These arguments render the 
validity of a full microscopic approach doubtful. Moreover, microscopic based relations 
can be relatively complex and use a significant number of parameters to describe the 
material involved. The main goal of this article is to find an improved scaling relation 
which balances the need for better founded formalisms and practical applicability to 
measured results in wires. As a starting point we use three central statements, which 
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will be detailed further throughout the article. 

Our first statement is that the temperature dependencies Hc2{T) and H^iT) (the 
thermodynamic critical field) are known. We showed earlier [56,60] that the shape of 
Hc2(T) is constant for all available NbsSn results, independent of the sample layout, 
morphology, A15 composition, strain state or applied critical state criterion and also 
holds for inhomogeneity averaged, extrapolated, H*2{T) data. This means that the 
normalized temperature dependence of Hc2 can be described by a single, known function 
from the microscopic theory. Also Hc{T) is known [11,68,69] and k,i{T) immediately 
follows through ki = Hc2/ {\/'2Hc) [29]. This approach differs from the temperature 
dependence in the Summers relation, in which Hc2{T) is calculated from an estimated 
K,{T) [1]. 

Our second statement is that non-linearities, which generally appear in Kramer 
plots [16] when approaching the maximum Hc2 that is present in a wire, can be attributed 
to A15 inhomogeneities. This assumption can be validated for specific wire layouts 
through a combination of simulations and measurements [60,62]. Alternatively, one 
can chose to fit measured non-linearities assuming a different pinning model as is often 
done [5,7,15,70]. This will improve the quality of the fit per specific wire type when 
approaching the maximum Hc2, but this is also the regime where Jc becomes very low. 
This regime is therefore of limited interest for applications. Fitting this regime per wire 
type will obviously result in a reduced overall error but also in a model with significantly 
reduced generality. We therefore chose not to include such non-linearities in the model 
fit to arrive at a generally valid relation for all NbaSn wires, with a validity up to the 
regime where the effects of A15 inhomogeneity start to appear (typically from about 
80% of the maximum Hc2)- 

Our third statement is that the introduction of strain sensitivity occurs solely 
through strain induced changes in Hc2, and ni. Alternative approaches can lead 
to a strain dependent pre-constant in the Jc{H,T,e) relation [5-7,70], which implies 
that the pinning efficiency changes with strain. Although this cannot be ruled out, it 
can be argued that the required change in pinning efficiency of at least a factor of 2 
over the relevant strain regime (below 1% axial strain) [7,70], is unrealistically large. A 
more reasonable approach to our opinion, is to assume that the main effect of strain is a 
modification of the electron-phonon interaction spectrum, which is supported by recent 
calculations [11,12,24-26]. This will obviously change directly T^, Hc2 and He and thus 
Ki. We will show below how this assumption, combined with a known temperature 
dependence of Hc2 and leads to a consistent introduction of strain and temperature 
effects in the Jc{H,T,e) relation by stating that H*2 — ^ H*2{T,e) and ki — > ki(T, e). 
Note that this is contrary to the usual unification of strain and temperature dependence, 
which assumes two separate and different dependencies on H*2(T) and H*2{e) [71]. 

We will show below how these three starting points lead to a consistent, simple 
and accurate relation for Jc{H,T,e). In Section |21 we will summarize the descriptions 
that are available throughout the literature and arrive at a general form for Jc{H, T, e). 
In Section El we will systematically compare this general form to measurements to 
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determine the resulting scaling relation. The overall accuracy will be demonstrated in 
Section ^ and discussed in Section |S1 together with a determination of the minimal 
required set of measurements that are needed to fully characterize a specific wire. Our 
conclusions are presented in Section IHl 



2. The critical surface 
2.1. Magnetic field dependence 

2.1.1. Relations for the Pinning Force The magnetic field dependence of the critical 
current density in NbaSn is determined by de-pinning of the fiux-line lattice and thus 
by the magnetic field dependence of the bulk pinning force Fp = ~Jc x B. The latter 
depends on the microscopic de- pinning mechanism. A general description for Fp{H) 
results from the observation that, for a multitude of low temperature superconductors, 
the bulk pinning force scales according to: 

FpOC^/(/i), (1) 

as was first pointed out by Fietz and Webb [72]. In (jT)), /i is the reduced magnetic field 
H/Hc2- The value for u is often obtained by neglecting ki and plotting the maximum 
bulk pinning force against Hc2 on a double logarithmic scale. Values for i/ found in 
this way are usually between 2 and 3 and for NbaSn values of 2 and 2.5 are mostly 
cited. However, since ni may change up to 50% with temperature f Section |2.2.3|) . its 
temperature dependence cannot be neglected for an accurate description of the critical 
surface of NbsSn. 

Several models for f{h) exist [16-18], but most often the Kramer version is used. 
Kramer derived two models for f{h), using two regimes. For the regime below the peak 
in f{h) he assumed that flux motion primarily occurs by de-pinning of individual fiux- 
lines, whereas in the regime above the peak in f{h), de-pinning occurs by synchronous 
shear of the fiux-line lattice around line pins which are too strong to be broken. The 
latter is referred to as the flux shear model and for this regime Kramer proposed [16]: 

Fp = ^ (2) 

127r2 (l - ao^/p) flo 

where Cee represents the elastic shear stiffness of the fiux-line lattice, oq is the fiux-line 
spacing and ^/p represents the density of the pinning planes. In bulk NbsSn, the grain 
boundaries are the main pinning centers [73-77] and can thus be interpreted as the 
pinning planes yielding l/^/p = c^av, the average grain size. 

One relation for the shear modulus of the fiux-line lattice of a high field 
superconductor without paramagnetic limiting, was derived by Labusch [20]: 



Cee{H) = 7.4xlO'^l^^^^{l-hf [N/m^] . (3) 



tc2l 

K 



"1 
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The flux-line spacing can be estimated by assuming a triangular flux-line lattice [68]: 

//, \ 0.25 / , \ 0.5 

-<-'^a) (^) • « 

Combining Q with (jH)), assuming l/i/p = d^v and a triangular flux-line lattice with a 
spacing as given by (jH) results in: 

Fp = 12.8 ^^°^f^''' ^° '(^-^)' ^ [GN/m=^] . (5) 



1 1^1 - Oa (i^)/(iav 



Using the Kramer form therefore results in u = 2.5 and 7 = 2 in (Q), if /(/i) is defined 
as f{h) = h'^-^{l — hY and (1 — a ^{H) / dg_-y) = 1, which will be validated below. Many 
authors indeed state u = 2.5 [1,16,17,62,71,78]. 

Higher values for u can be found if experimental results are described with non- 
Kramer-like pinning. Deviating values for i/ in are reported by Kroeger et al. [70] 
who claim i> = 3.0 and by Cheggour and Hampshire [5] who state u = 3.2. It should 
be noted that both studies neglect ki and fit the experimental results to Fp oc (1 — /i)^ 
and Fp oc (1 — h)^'^ respectively. This choice suggests different pinning behavior. It will 
be explained below that such high values for the power above {1 — h) can arise when 
sample inhomogeneities are ignored. Here, they will be accounted for, which results in 
the experimentally supported assumption that Fp oc (1 — /i)^, as will be discussed in 
more detail in Section l3?Tl 

Combining ^ with the reverse definition of the Lorentz force gives (Suenaga and 
Welch in [79]): 

x0.5/„ o-^0.25 1-1x10' /Xo(gc2-g) 

Ki (1 - (f/j/dav) 
In most technical wires d^v is 100 to 200 nm, whereas the flux-line spacing decreases 
rapidly below 50 nm for magnetic fields above 1 T. Therefore ai^{H)/da_w -C 1 and 
Jc'^if^oH)^'^^ will approximately be linear in H and can be referred to as a Kramer 
function /k: 

1 0^ 

/k (H) ^ J°-^ if^oHf' - —^0 - H) . (7) 

A linear extrapolation of ((Tj) to Jc = then yields Hc2 and the slope can be used 
to calculate ki. The Kramer extrapolated critical field will be referred to as to 
distinguish it from a measured Hc2- 

The pinning force, as given by peaks at /i = 0.2 if a^{H) <^ rfav to yield the 
maximum pinning force: 

Fpmax = ^-^^^f"^' ' [GN/m^] . (8) 
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2.1.2. Low field linearity of Kramer plots The dependence Fp oc (1 — hY occurs when 
shear deformation of the flux-hne lattice is the primary de-pinning mechanism as was 
mentioned above. This was shown experimentally by Cooley et al. [80] to be the case 
when the grain size is substantially larger than the flux-line spacing (c/av > 2a/x(-ff)). If 
the grain size is comparable to the fiux-line spacing, a direct summation of individual 
pinning interactions between grain boundaries and fiux-lines occurs, which corresponds 
to Fp oc (l — h) [20,22,80]. Plots of the Kramer function /k from ((7j) versus magnetic field 
(a 'Kramer plot'), calculated from Jc data measured on NbsSn conductors, are generally 
linear in the 'standard' Jc measurement regime of about 0.2 < h < 0.6. This supports 
shearing of the fiux-line lattice as the primary de-pinning mechanism, provided that 
the grain size is larger than about 100 nm [80]. Moreover, Kramer plots resulting from 
magnetization data on Powder-in- Tube processed wires, reacted at various temperature 
and time combinations, were shown to be perfectly straight down to h = 0.07 at 
various temperatures by Fischer [58]. Also transport Jc measurements on different 
bronze processed wires indicated straight Kramer plots at various temperatures down 
to h = 0.04 [10,78]. Observed downward curvature when approaching lower magnetic 
fields can, apart from a different pinning behavior, additionally arise through self-fields 
generated by the large transport current densities. Accurate correction for these self- 
fields is non-trivial. 

These observations validate the application of © to very low fields and support the 
assumption that the (1 — a^/fiav) term has negligible infiuence and can be approximated 
by (1 — a^/iiav) — 1- Moreover, they indicate that for technical wires, which have an 
average grain size above 100 nm, shear deformation of the fiux-line lattice apparently 
remains the primary de-pinning mechanism also below the maximum in Fp. The 
assumptions used to arrive at the Kramer model have received some criticism [21]. 
In practice, however, using the high magnetic field description of Kramer appears to 
hold for NbaSn wires, at least down to 1 T. 

2.1.3. Linearity of Kramer plots approaching Hc2 Above h = 0.8, non-linearities 
are mostly observed in Kramer plots of wire results when approaching the maximum 
Hc2 that is present in the wire. The pinning force function combined with the 
approximation (1 — a i^{H) / d^^^) = 1 is therefore often written in a more general form: 

Fp (H) = Fp^^Jih) = C^J^^^h^ (1 - hy , (9) 

where C is a constant and p and q determine f{h). This allows for insertion of different 
values of p and q that implicate a different pinning behavior, but most probably originate 
from inhomogeneity averaging. Such different values for p and q lead to still ongoing 
discussion in the literature whether or not a general pinning mechanism, or f{h), can 
be assumed for NbsSn wires. 

The field dependence of the bulk pinning force as given by (jH)) is shown in figure C] 
and combined with Kramer plots for various values of q around 2. A variation of p only 
infiuences the low field region. The aforementioned observations of perfectly straight 
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Figure 1. Calculated changes with q in the field dependence of the bulk pinning force 
(left plot) and Kramer function (right plot). 

Kramer plots at lower magnetic fields for bronze and Powder-in- Tube processed wires, 
however, render strong deviations of p from 0.5 unlikely. 

A ± 25% variation in q results in a change in the maximum pinning force and 
its position and an up- or downward curvature in the Kramer plot for h —>■ 1. The 
effect of a lowering of q, i.e. an increase of the maximum pinning force combined 
with a shift of its position towards higher magnetic field is identical to what has been 
observed experimentally through grain refinement in Powder-in- Tube wires by Cooley 
et al. [80]. Upward and downward tails in Kramer plots have both been observed 
experimentally (e.g. Ekin in [15], [60], Suenaga in [81], [82]) and have been attributed 
to microstructural [83] and compositional (Ekin in [15], [60], Cooley et al. [62], Suenaga 
and Welch in [79]) origins. A range of experimentally detectable Hc2 values could 
also be attributed to flux-creep and thermal activated fiux-fiow (e.g. the presence of an 
'irreversibility field' in high temperature superconductors), but we find that these effects 
are negligibly small compared to the severe influence of compositional inhomogeneities 
on Hc2 [64], especially at temperatures below 20 K. 

The most convincing support for attributing curvature in Kramer plots when h ^ 1 
to inhomogeneities is given by Cooley et al. [62] , who numerically modeled the filaments 
in a Powder-in- Tube processed wire as concentric shells with different Sn content and 
thus different Hc2 and Tc. It was assumed that each composition behaves according to 
q = 2. An area weighted summation of the transport current density of such a system 
of parallel paths of variable property showed that a concave tail appears in the overall 
Kramer plots when the maximum available Hc2 is approached. Increasing the severity 
of the Sn gradient in the A15 sections increases the occurrence of such tails in otherwise 
linear overall Kramer plots. This parallel path model is supported by measurements 
on similar PIT processed wires [56,60]. Such an area weighted summation over parallel 
paths of lesser and better A15 quality results in a similar overall Kramer plot as for 
p = 0.5 and q = 2.5 in figure ^ A downward tail in Kramer plots, as is often observed 
in bronze processed wires, can be described with g < 2 as is shown in figure ^ thereby 
again implicitly assuming a different de-pinning mechanism. However, assuming a serial 
connection of lesser and better A15 quality regions for the more irregular filaments in 
bronze processed wires, also results in a downward tail in Kramer plots when /i — 1, 
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since the lower quality A15 sections will determine the overall transport properties. This 
argumentation is a strong indication that the curvature is not due to different de-pinning 
mechanisms but to inhomogeneity effects and thus that q should be 2. 

It can thus be stated that the Kramer ffux-line shear model, at constant 
composition, is valid over nearly the full magnetic field range, i.e. up to Hc2 and 
down to a magnetic field where the flux-line lattice spacing becomes approximately 
more than half of the A15 grain size [80], which is below 1 T from (0)). Non-linearities 
in Kramer plots when approaching the maximum present Hc2 are thus attributed to 
inhomogeneity averaging. The effective bulk values H*2 and T*, as detected by critical 
current measurements, are weighted averages over the A15 volume, but it is not always 
clear how the averaging occurs. The discussion above leads to the conclusion that there 
is no a priori reason to deviate from Fp oc h^'^{l — h)'^, provided that the inevitable 
presence of A15 inhomogeneities in wires is recognized. 



2.2. Temperature dependence 

To come to an overall description for Jc{H, T, e), the temperature dependence of the bulk 
pinning force is considered. In observing Q it is clear that temperature dependence 
occurs through temperature dependence of Hc2 and ki. 

FpiH,T)^C^-l^^^^^^fih), (10) 

with h = H/Hc2(T). It is thus required to discuss the available descriptions for Hc2(T) 
and Ki{T). For both parameters, empirical as well as microscopic alternatives are 
available. In the Summers relation [1] empirical forms are used. It will be shown here, 
how these can be replaced by alternatives with a better connection to the microscopic 
theory. 



2.2.1. Empirical temperature dependence First, the empirical temperature dependence 
of the upper critical field, as used in the Summers relation, is introduced. Based on 
earlier work by Suenaga in [81] and Hampshire et ai. [78], Summers et ai. [1] proposed 
an empirical relation for the temperature dependence of the upper critical field, thereby 
improving the correspondence between the available scaling relations and measured 
results: 

Hc2 iT) _ ^2^^liT) 
H,2 (0) ^ ' (0) 

= (1 - t^) [1 - 0.31t^ (1 - 1.771nt)] , (11) 

where t = T/T^. In this relation an empirical fit for Ki(T)//ti(0) is inserted in 
Ki = Hc2/{y/'2Hc) and Hdt) / Hc{0) = (1 — t^) is used for the temperature dependence 
of the thermodynamic critical field. It has to be emphasized that Hampshire et ai. [78] 
claim a validity range only below 13.5 K for their calculated ki(T) values. 
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2.2.2. Microscopic based Hc2{T) An early microscopic description for the temperature 
dependence of the upper critical field was simultaneously derived in 1964 by Maki [28,29] 
and De Gennes [27,30] which, written in the De Gennes form yields the implicit relation 
(MDG relation): 

The function uses only two parameters, namely Tc(0) (which can be measured) and 
the diffusion constant of the normal conducting electrons D. The diffusion constant is 
inversely proportional to the slope of the Hc2(T) dependence at Tc{0) [56,60]. The other 
parameters are the reduced Planck constant (h), the magnetic flux quantum (0o) and 
the Boltzmann constant (/cb)- The terms ^p{x) represent the digamma function [84]. 
The description was derived assuming a dirty superconductor (i.e. i -C it uses a 
spherical Fermi surface approximation and assumes a constant density of states at the 
Fermi level N[Ep) and a weak electron-phonon interaction. Furthermore it assumes no 
paramagnetic limitation of Hc2{T) and an absence of spin-orbit scattering. A normalized 
form of ()12|1 with respect to temperature can be defined: 

M/^G(t) = §^i^, (13) 

in which i^c2(i':)MDG represents Hc2(T) calculated using (|T^ . 

It can be convenient for practical applications to avoid the need to solve the implicit 
MDG relation by replacing it with an explicit expression. An approximation of (fT^ for 
the entire temperature range yields: 

MDG (t) = (1 - t^-^^) . (14) 

This power of 1.52 is practically identical to the value of 1.5 which was found 
from a similar fit as (fT^ to (partly) extrapolated H*2(T) results by Cheggour and 
Hampshire [5]. A power of 1.5 was also found recently from analysis using Eliashberg 
theory by Oh and Kim [11]. 

Application of the full Eliashberg based formalisms is the most correct approach, 
since it accounts for all the electron-phonon interaction strengths. However, there are 
two reasons for applying only the simplest form of the microscopic descriptions [i.e. 
(I12|l ] to wire results. First, since wires are inherently inhomogeneous there is no single 
Hc2{T) but a distribution of properties [60]. Also, the inaccessibility of the normal state 
resistivity (which is a prime parameter in the microscopic theory) of the A15 in wires 
renders a proper connection to the theory practically impossible. Second, the use of 
the formal descriptions results in an increased number of parameters compared to (fT^ . 
This can be avoided by implementing the additional parameters that are required for a 
proper connection to the microscopic theory in D [and thus in the slope at Tc(0) [60]] 
and Tc(0) in (fT^ . To allow this it has to be assumed that the shape of Hc2{T), at least 
to first order, remains unchanged. 

The simplest form of the microscopic descriptions, the MDG relation (fT^ . was 
tested previously and found to fit Hc2{T) for all available literature NbaSn results and 
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measurements on a multitude of wires with very high accuracy [60]. This accuracy 
holds for the Orlando et al. thin film results [48] at all Sn concentrations, the Foner and 
McNiff single- and polycrystalline results [50,51], and the Jewell et al. bulk results [67]. 
The single crystals are close to stoichiometric and thus strong coupling [63] and most 
probably clean. The thin films and bulk samples are of varying degree of resistivity 
and Sn content and cover therefore a large range from approximately clean to dirty 
and weak to strong coupling. This leads to the conclusion that the shape of the field- 
temperature phase boundary is, within the experimental error bars, independent of the 
electron-phonon coupling strength (which can also be confirmed from the theory [40]), 
independent on whether the material is clean or dirty, and not influenced by the 
(unlikely [44, 48, 56, 85]) presence of PPL and spin-orbit scattering. It can thus be 
expected that it can be validated for the entire range of compositions that are present 
in wires. This is confirmed by the fact that it also holds for Kramer extrapolated, 
inhomogeneity averaged H^(T) results [60]. 



2.2.3. Microscopic based ki(T) From the microscopic theory [41] it follows that for 
weak coupling superconductors in the dirty limit Ki(0)/ki(Tc) = 1.2. In the strong 
coupling limit this ratio saturates at about 1.5 [41], which means that the maximum 
expectable change in ki(T) is 50% from K to Tc. 

For calculation of ki(T)/ki(0) for a given interaction strength the nonlinear 
Eliashberg equations [13] will have to be solved as was done numerically by Rainer 
and Bergmann [41]. No general simple function for /ti(T)//ti(0), derived directly from 
the electron-phonon spectrum and thus valid for all Aep, can therefore be given for 
wires. Fortunately, a generalized function can be derived for fi;i(T)/Ki(0) through 
Ki = Hc2/{\/'2Hc), using the microscopic form for Hc2{T) / Hc2{0) (fT^ combined with 
the temperature dependence of He- 

It is well accepted that Hc{t)/Hc{0) is for all superconductors very close to 
(1 — t^) [68]. Recent calculations based on Eliashberg theory yield Hdt) / Hc{0) = 
(1 — t^'^^) [11]. The deviation from (1 — t^) is measured for NbaSn [69] and is within 
about 2%. A fit to the measured Hc(t) / Hc{0) yields (1— 1^-°^). For practical applications 
a power of 2 will thus be sufficiently accurate. Combining (fT^ / p3|l or p4|l with 
= H^^/iV^H^) and H^{t)/H^{0) = (1 - t^) thus yields: 

Mtl^MDGW^l^, (15) 

^ ^ (0) 1 - 12 1-^2 V ; 

By implementing an accurate microscopic form for Hc2{T) in the scaling relations also 
Ki{T) is accurately known, rendering an empirical version redundant. 



2.3. Strain dependence 

2.3.1. Strain dependence in Nb^Sn wires When a NbaSn composite wire is subjected 
to longitudinal compressive strain, the critical current density, upper critical field and 
critical temperature reduce approximately linearly and reversibly with strain. When 
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a wire is subjected to a longitudinal tensile strain, the critical parameters increase 
approximately proportionally and reversibly with strain until a parabolic-like peak is 
reached, after which the properties reduce approximately linearly with strain. The 
larger thermal contraction of the matrix materials compared to the thermal contraction 
of the NbaSn results in an axial pre-compression of the A15 when cooled from reaction 
temperature to the test temperature below 20 K. When the wire is axially loaded after 
cool-down, this pre-compression is minimized and J^, i?c2 and increase. The position 
of the maximum in the strain dependency curve appears at the point where the three- 
dimensional deviatoric strain components in the A15 are minimal. In practice this 
correlates closely to the point where the axial pre-compression of the A15 is minimal. 
The axial pre-compression is mostly identified by em, whereas the position of the 
minimum in the deviatoric strain components is identified as 6. Several models have 
been introduced in the literature to describe this behavior, ranging from polynomial fits 
to relations based on the Eliashberg theory. These will be discussed below. 

When an A15 lattice is deformed, its vibration modes will change. In addition, 
the electronic structure will be modified and therefore N{Ef). It is thus reasonable 
to expect a strain induced variation of the electron-phonon interaction spectrum and 
the phonon density of states and hence, a change in the interaction constant Aep. The 
difference between the strain sensitivity of Hc2 and represents an important link to 
microscopic understanding since both depend differently on N{E-p) and Aep. It is found 
experimentally that [71,94]: 



where the index 'm' indicates the vahics of Hc2 and at the maximum of the axial 
strain dependency curve. The power zu = 3, indicates that Hc2 is roughly 3 times more 
sensitive to strain than Tp. This was found to be vahd at all temperatures (e.g. [3,4]). 
An attempt was made by Welch [94], following earher work by Testardi [95-97], to 
explain the observed difference in strain sensitivity on the basis of the strong coupling 
renormahzed BCS theory, using a McMillan [98] based, Allan and Dynes [99] strong 
coupling formulation for Tc. It was concluded that w = 3 cannot be explained by 
either a change in N{E-p) or Agp alone and the strain dependence should therefore be 
described by a combination of both. It should be pointed out, however, that McMillan 
based descriptions for Tc are valid for Aep < 1.5 [100], whereas reported values for 
Nb3Sn are mostly higher [48] . It is clear, however, that improved understanding of the 
value for zu should result from exact calculations of strain induced modifications on 
the full electron-phonon interaction spectrum and the density of states. This might be 
possible through a combination of the recent efforts of Markiewicz [24-26] and Oh and 
Kim [11, 12]. An important conchision drawn by Welch is that in a three-dimensional 
strain description the deviatoric strain components dominate the strain sensitivity of 
NbsSn and that hydrostatic components have, in comparison, a negligible effect. 

Most strain dependencies are defined through the strain dependence of the upper 
critical field or the critical temperature. Two general strain dependent terms will 



H,2 (4.2ir,e) 

iJe2m (4.2ir) 




(16) 
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therefore be defined: 




(17) 



under the assumption that s(e) is independent of temperature and zu ^ 3. 

2.3.2. Available models Five models describing the strain sensitivity of NbsSn can 
be distinguished throughout the literature. The simplest form is an introduction 
of hydrostatic strain though a pressure term. Although this yields a relatively 
simple connection to thermodynamic calculations to describe single crystal hydrostatic 
experiments, it is not suited for more practical systems such as wires due to the 
dominating deviatoric components. An empirical fit for axial deformations in wires 
was recently adapted by Hampshire and co-authors [6] in the form of a fourth-order 
polynomial. Although this will obviously result in accurate fits to experimental data, it 
lacks any connection to the underlying physics, and is thus not suited for the goals set 
out for this article. The three remaining models that have been proposed to describe 
experiments on wires will be discussed next in chronological order. 

The first model that was proposed to describe the strain dependence of the critical 
properties of NbsSn composite wires was the so-called power-law model, introduced by 
Ekin in 1980 [71]. The model was developed based on measurements using an axial 
pull strain device. The investigated compressive strain range was therefore limited to 
the thermal pre- compression of the A15, introduced by the matrix components of the 
wires. A power-law dependence was proposed for s(e) using the extrapolated upper 
critical field [71]. In this power-law dependence, a different axial strain (ca) sensitivity 
is used for ea < and for eg, > 0. The power-law model is an empirical description 
that does not account for the three-dimensional nature of strain, but includes the 
generally observed asymmetry in axial strain experiments. Observable linearities at 
large compressive longitudinal strains are not accounted for and it is therefore valid 
over a limited strain regime. 

The second model that was proposed was the so-called deviatoric strain model. For 
high compressive strains Ten Haken found an approximate linear dependence of on 
axial strain, deviating from the power-law behavior [102, 103]. These deviations from 
the power-law model at large compressive axial strain values, combined with results of 
experiments on quasi two-dimensional NbaSn tape conductors [104-106] (which allow an 
analytical calculation of the three-dimensional strain components) led to an alternative 
description based on the non-hydrostatic, i.e. distortional or deviatoric strain [8,104]. In 
this deviatoric strain model only the second strain invariant is considered. The model 
is empirical, does account for the three-dimensional nature of strain and linearity at 
large compressive axial strains, but is symmetrical around the maximum in the strain 
dependency curve, and is claimed to be valid only for the compressive strain regime. 

The third model available in the literature is the so-called full invariant strain 
analysis [24-26]. This model is based on the assumption that the strain dependence 
of Tc results from the strain dependence of the phonon modes in the NbsSn crystal 
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lattice. The model uses a full invariant strain analysis which correlates the strain 
dependence of Tc to strain induced changes in the phonon frequency spectrum using 
relations that stem from Eliashberg theory and a strain energy potential function. The 
main result of the full invariant model is that it is able to demonstrate the effect of the 
harmonic and anharmonic terms in the strain energy potential function independently. 
It explains how the reduction due to hydrostatic strain decreases linearly with axial 
strain. It further explains the experimentally observed parabolic-like behavior around 
the maximum, the linearity for large compressive axial stain values and the asymmetry, 
in terms of the second and third strain invariants. For now, unfortunately, it does not 
include a description for H*2{e), which is required for Jc scaling. The latter might be 
provided by a combination of recent promising work by Oh and Kim [11, 12] and a full 
invariant analysis. 

2.3.3. Empirical correction to the deviatoric strain model The full invariant analysis is 
the most proper method available for the introduction of strain dependency of the critical 
properties. The resulting relations, however, remain somewhat complex compared to 
the empirical fits, and a complete description for Tc(e) and H*2{e) is not fully developed 
yet. It is clear that the empirical models all lack specific details in comparison to the full 
invariant analysis. An empirical model that accounts for three-dimensionality, linearity 
at large axial strain, and asymmetry, in agreement with the invariant analysis, can be 
obtained by assuming a simple mathematical correction to the Ten Haken deviatoric 
strain model [56]. Our starting point is relation (3) in [8]: 

f^oH^ (cdcv) = l^oH^ (0) — Cdcv\J (cdcv)^ + (eo,d)^, (18) 

in which edev represents the second deviatoric strain invariant, Cdev represents the slope 
in the hnear regime of /io-f^^(edev), and eo,d represents the non-axial, remaining strain 
components when the axial strain is zero. In axial form (fTHj) can be written as [8,56]: 

^i,Hl (6a) = ^l,Hl (0) - C[^J {e^f + {eo,.)\ (19) 

where C'^ is a constant and eo,a represents the remaining strain components, as with the 
factor eo,d in (fTHj) . To account, in the axial form of the Ten Haken model, for a stronger 
reduction at tensile strains, a linear term causing a reduction of with axial strain is 
required. This can be achieved by including an additional linear term in p8|) : 

/^o-f^,^ (cdov) — (0) — Ci^^2\j (cdov)^ + (eo,d)^ 

~ CinvS^dev, (20) 

where Cinv2 and Cinvs are constants representing the second and third strain invariants. 
The term eo,d now represents the initial suppression of due to hydrostatic strain. A 
similar route can be followed for the axial version ()19p . resulting in: 

li,Hl (ea) = ^^oHl (0) - C'^.^J {e^f + {e.^.f - C',,e.,. (21) 



A general scaling relation for Nh^Sn 



14 



in which ea = eappUod + ^- The correction C^2^a effectively rotates the strain 
dependency curve around its maximum. This causes the axial position of the maximum 
to shift by a factor esh = Cl2^o,a/[(Cli)^ — {C'^2)'^T'^ ■ large asymmetry, 5 (the 

axial strain at which the minimum in the deviatoric strain occurs) will then deviate 
significantly from the observable axial position of the maximum em- In axial strain 
experiments it is more convenient to normalize the strain dependence function to the 
observable maximum, also to avoid the introduction of an additional parameter through 
the use of 5 in combination with Cm- The rotation around the maximum also causes 
the normalized form ()17|) to increase above 1 at the position of the maximum. This is 
undesirable since in scaling relations the reduction in e.g. is usually normalized to 
the maximum observable value at the peak of the measured strain dependency curve, 
and not to the 'strain free' value since this is not accessible in an axial strain experiment. 
This can be prevented by re-normalizing s(e) to the maximum observable critical field 
value. A re- normalization of s(e) to the observable peak position and peak value in an 
axial strain experiment leads to a new axial form: 



It has to be emphasized that this is an empirical correction to include the effects of 
the third strain invariant which, as suggested by the full invariant analysis, should be 
accounted for to describe asymmetry. It results in a simple description with three 
parameters (plus pre-strain) that accounts for the overall behavior as expected on 
the basis of a full invariant analysis and experimental observations. The hydrostatic 
strain components are included in the remaining strain term eo,a and the strain induced 
behavior resulting from the second and third invariants is accounted for through Cai 
and Ca2 respectively. Relation will be used through the remainder of this article 
since it is the simplest form that still accounts for the three-dimensional nature of strain 
and is based on an underlying physics model. Possible future, more fundamental based, 
strain descriptions can be easily implemented in the general scaling relation, since all 
strain dependence is included explicitly in s(e) in terms of a change in Hc2 and/or 





(22) 



2.4- Field, temperature and strain dependence of the bulk pinning force 



The strain dependence described above has to be included in the relation for the 
magnetic field and temperature dependent bulk pinning force pO|) to arrive at a relation 
that can be used to inter- and extrapolate critical current measurements on wires. This 
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requires the use of a strain dependent critical temperature as defined earlier by (fTTj) and 
a temperature and strain dependent GL parameter and upper critical field: 

H,2{T,e)=H,2m{T)s{e). (23) 

Implementing this in the Kramer form of the field dependence of the bulk pinning force, 
i.e. (0), and assuming (1 — a^/c/av) = 1 yields: 

Fp {H, T, e) = ^^'f (1 - hf . (24) 

Ki (T, e) 

The upper critical field thus delivers a term s(e)^'^ acting directly on the bulk pinning 
force in addition to strain dependence arising through f{h) and Tc(e). The strain 
dependence of ki is, however, undefined since this requires knowledge of the strain 
dependence of He which can differ from ifc2(e)- In addition, some strain dependence 
through strain induced changes of the fiux-line to lattice interactions could be assumed. 
To arrive at an overall description, it is thus required to have knowledge of ki(T, e), 
and/or to introduce a strain dependent constant C(e). 

2.4-1 ■ Ekin's unification of strain and temperature dependence In analyzing critical 
current versus strain data at 4.2 K Ekin [71] found that the bulk pinning force scales, 
in analogy to the Fietz and Webb scaling law (^, as: 

F^{H,e),,^o:[H:,{e)ff{h) (25) 

where n = 1.0 ± 0.3 and the star again indicates a bulk average. The power of 1 differs 
from the power of 2.5 that is required for magnetic field and temperature scaling, as 
seen from (jlOp . Ekin solved this inconsistency by postulating a combination in terms 
of explicit dependencies on temperature and strain. Neglecting the temperature and 
strain dependence of ni he proposed: 



F^{H,T,e) =Cs{e 



f{h) 



L^c*2(0,e) 

^Cs{er{l-tYf{h) (26) 

where h = H/H*2{T,e), t = T/T*{e) (stars indicating bulk average values), n = 1, 
V = 2.5. Note that the third term only describes the normalized temperature dependence 
of the upper critical field approximated as (1 — t^), in which strain sensitivity occurs only 
through T*{e). The combination of strain and temperature sensitivity was introduced 
on a fully empirical basis, based on 4.2 K critical current data. This implicates, as 
was pointed out by Ekin, that the term C7//«i(T, e)^ in ()24|) has to produce a strain 
dependency s(e)^'^ to counteract the s(e)^'^ term arising through the strain dependent 
upper critical field in ()24|1 . A similar magnetic field, temperature and strain dependence 
is used in various forms of generalized scaling relations which, by combining ()26p with 
Fp = — J X B, result in an overall description for the critical current density in NbaSn 
wires. 
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2.4-2. Unification of the Kramer form The ratio [fioHc2{T,e)]'^'^/Ki{T,ey in (j24p can 
be normalized to zero temperature to express the direct strain influence on Fp (i.e. not 
through Tc(e)) separately, so that the result only contains temperature dependence and 
strain dependence through Tc(e). With the definitions: 

this leads to: 

T? (TT T A - [^0^^2m (0)]^'^ ^ ^^^2.5 

, ej ——2 s {e) 

Ki (0, e) 

X HT,nie)f' _ ^)2 _ (29) 

To retain consistency with Ekin's observation that -Fp(e)4.2 k oc s(e) (under the 
assumption that this is valid for all temperatures) it can for example be postulated 
that: 

/^i(0,e) = /ci^(0)s(e)^, (30) 

leading to: 

i2.5 



F,(g.r.e)= ^l^°^';°'°>l' .(e) 

Ki (0) 



^ bjT^TUe)r^o., . _ (31) 

kiT,T,ie)f 

Alternatively, one can define a strain dependent pre-constant C{e) as will be shown 
below. The temperature dependencies in fl27|) can be empirical (e.g. (fTT|) ) or microscopic 
based alternatives (e.g. (fT^ and (fT5|) ). The discussion above is used in the scaling 
relations that are available throughout the literature. 

2.4.3. The Summers relation Summers et al. [1], following previous work by 
Hampshire et al. [78] and Ekin [71], using e = e^, defined a strain dependent pre- 
constant: 

C{e) = Cos{ef'. (32) 

The power of 0.5 is introduced in ()32j) to retain consistency with the strain dependence 
of the bulk pinning force as postulated by Ekin ^I^i . It can be shown that the Summers 
relation can be rewritten in a form consistent with the previous Section as: 

Fp (ff, T, £) = C„ |fi„ft2„ (0)1° ' s (e) 

k{T,T^ (e)) 

in which the empirical temperature dependence (fTT|) is used for b(T,Tc{e)) and 

k{T,Tc{e)). The Summers relation is thus, apart from the pre-constant, consistent 
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with (|^. It has finally to be quoted that Hampshire et al. stated [78], in relation to 
substantial inhomogeneity of their samples, that the validity of their calculated /«i(T) 
results is expected to break down for T > 13.5 K. It will be shown below that this, 
through comparison of the Summers form (fTT|) to the microscopic based form (fT3j) and 
measured results, indeed appears to be the case. 

2.4-4- Alternative approaches The empirical bases in the relations proposed by Ekin 
and Summers et al, their inaccuracy, and the extensive new Jc{H,T,e) data sets for 
wires that have become available over the past decade have evoked a re-analysis of overall 
scaling behavior in NbaSn wires. Two alternative approaches can be distinguished, 
focusing either on improving the strain and temperature dependent terms [4,9], or 
focusing on improving the basis of the strain dependence of C//ti(T, e)^ in [5-7]. 

In the first approach it is stated that, since the strain dependence through the 
power-law function does not account for linearity at high compressive strains, it should 
be replaced by a function that does. This was done through a change in s(e) from 
the power- law form to the Ten Haken form. Secondly, it was assumed that (j3T|l . and 
thereby the Summers version, is valid but that the empirical temperature dependencies 
as proposed by Summers et al. [i.e. (fTT|) ] are in error and can be improved upon. 
The general form ()1()|) was chosen in [4,9] under the condition that s(e) from ()17|1 is 
independent of temperature, (fT^ has general validity and (1 — a/^(i/)/(iav) — 1- 

The temperature dependence in (fTIH) is determined by the ratio [fioHc2iT)Y / niiT)'^ . 
Four- dimensional parameter least squares fits using extensive data sets on a six different 
wires and highly linear Kramer plots indicated n = 1 (as stated by Ekin), p = 0.5 and 
g = 2 (as stated by Kramer for h > 0.2). Since the temperature dependencies of Hc2 and 
Ki were fully empirical in the Summers form, both powers u and 7 were regarded as free 
parameters to allow for errors in the temperature dependencies. Overall least squares 
fits on the data sets indicated that u = 2 and 7 = 1 yielded the highest accuracy. The 
overall description (fTO|l with p = 0.5, q = 2, u = 2 and 7 = 1 yielded fits to experimental 
data with a standard deviation of about 2% over a limited magnetic field range from 5 
to 13 T, but for all temperatures (4.2 K to Tc) and all investigated compressive strains 
(—0.7 < ea < 0) [2]. It was expected that if correct dependencies for Hc2{T) and 
/€i(T) were used, that similar accurate fits would be obtained with values for p and 7, 
which are consistent with Kramer's pinning theory, i.e. u = 2.5 and 7 = 2. However, in 
Sectional we will re-analyze one of these earlier data-sets with the recent new insights on 
inhomogeneity averaging effects and the now known temperature dependencies Hc2(T) 
and Ki(T), and show that the deviating values u = 2 and 7 = 1 are correct. 

The second approach to improve upon the Summers relation emphasizes on the 
different dependencies of the bulk pinning force on temperature and strain in an attempt 
to define a better founded combination of both dependencies. Following initial work by 
Kroeger et al. [70] this results in a general relation of the form [5-7]: 

Fp {H, T,e)=A (e) [fioH:, (T, e)]" (1 - hf , (34) 
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in which n, p and q are free parameters and A{e) is a function of strain alone. Several 
versions for A{e) have been proposed [also including ni(T,e)] in attempts to describe 
strain induced modifications of the flux-line to pinning center interactions and the strain 
dependence of ni. 

2.4-5. Selected general scaling relation In general it can be stated that much of the 
developments on overall scaling behavior remains empirical and it is unlikely that this 
empirical basis can be fully removed. As a final remark it should be mentioned that in 
a four- dimensional description with a large number of free parameters it often remains 
a somewhat subjective choice which parameters are fixed or allowed to vary, in addition 
to the choice of relationships that can be implemented. 

The overall scaling law that is proposed in this article will now be summarized. It 
is assumed that the critical current density scales according to plj) in the most general 
form: 

C [fioH:,iT, e)f 
fioH Ki (T, e) 

where h = H/H*2{T,e), in which H*2 represents the inhomogeneity (bulk) averaged 
critical field at which the critical current density extrapolates to zero. The microscopic 
based temperature dependencies are selected for the upper critical field (fTSj) and the GL 
parameter (fTH|) . This leads to the general normalized form: 



Je (H, T, e) = x7 '^^ (1 - > (35) 



1 (n T A ^ [A^0-^c2m (0)] 



X ^—^hP {1 - hf , (36) 



MDG {ty 
TW 

where t = T/T*{e). The term s{ey~°''^ arises from the strain dependence of the upper 
critical field through ()17p and an unknown strain dependence of the zero temperature 
GL parameter, defined as: 

ft:i(0,e) = ft:i^(0)s(e)°. (37) 

The form for s(e) is chosen to be the improved deviatoric strain model ()22j) . 

The temperature and strain dependent bulk average critical field in (jHBj) is given by a 
combination of the microscopic temperature dependence (fT^ and the strain dependence 
(El: 

(T, e) = H:^^ (0) MDG (t) s (e) , (38) 

and the strain dependence of the critical temperature is defined by (jl7|) . 

The constants p, q, z/, 7 and a will have to be determined through systematic 
comparison to experiments, and w = 3 is used in accordance with the experimental 
results of Ekin. Note that this approach significantly differs from the usual separation 
between the dependence of the bulk pinning force on H*2{T) and H*2{e). Here this 
is not assumed a priori and strain dependency is simply introduced through a strain 
dependent upper critical field and GL parameter. This relation will be compared to 
measurements in Section 01 to determine the constants and to test its validity. 



A general scaling relation for Nb3Sn 



19 



3. Comparison to measurements 

In this Section we will compare the general Jc{H,T,e) relation to earlier 

experimental data. Experimental details, as well as the existing data-set were described 
elsewhere [4,10,107]. The selected data-set for the renewed analysis are Jc{H,T,e) 
results which were measured on an International Thermonuclear Experimental Reactor 
(ITER) type, Furukawa bronze processed wire. A more detailed description of this 
wire is also given elsewhere [60]. The results of this wire are selected because of the 
demonstrated high reproducibility, and since it was analyzed extensively in various 
laboratories around the world. It should be emphasized that the results for this wire are 
not unique, but representative for the scaling behavior for at least five other wires from 
various manufacturers [4,9]. Jc{H, T) results were obtained on an ITER-type Ti-6A1-4V 
barrel at electric field criteria of Ec = 10~^ to 5 x 10""^ V/m. Jc{H,T,e) results were 
obtained on a Ti-6A1-4V U-shaped bending spring at = 5 x 10~^ V/m and on a 
Ti-6A1-4V circular bending spring ('Pacman') at Ec = 10~^ V/m. The ITER barrel 
and U-shaped bending spring samples were heat treated together, whereas the Pacman 
sample was heat treated separately. All data are corrected for resistive currents that 
run parallel to the superconductor at any voltage criterion and the self-field generated 
by the wire is also corrected for. These details can be found elsewhere [56]. 

First the field dependence will be analyzed in detail. This results in the values 
for p and q in flHUj) . Hl^iT) and F^^i^^{T). Then the temperature dependence will be 
analyzed in detail and it will be shown that the Kramer /Summers form fails to describe 
the temperature dependence by the use of z/ = 2.5 and 7 = 2 in (jSHl)- The use of these 
constants results in a temperature dependent variation of ni that is significantly larger 
than can be expected from the microscopic theory. Next the improved temperature 
dependence, through the use of different values for z/ and 7, will be introduced and 
validated. This is followed by an analysis of the strain dependence using the improved 
deviatoric strain model. By using correct values for v and 7, it is not required to 
separate the temperature and strain dependence of the bulk pinning force. Also the 
strain dependence of -Fpmax is determined from the measurements. Finally, the minimum 
required data-set to determine the critical surface of an unknown wire will be discussed. 

3. 1 . Normalized hulk pinning force 

To investigate the magnetic field dependence of the bulk pinning force, each combination 
of temperature and strain was attributed to a unique but unknown effective upper 
critical field H*2(T,e) and maximum bulk pinning force Fpms.x(T,e). The bulk pinning 
forces were calculated from the measured critical current densities through Fp{H, T, e) = 
\Jc{H, T, e)fiQH\. Kramer plots are linear for this wire as will be shown below, indicating 
q = 2 and p = 0.5 as was shown in Section 12.11 To analyze consistency with 
these values in (P) for all temperatures and strains, a least squares fit was made to 
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Figure 2. Reduced bulk pinning force as function of reduced magnetic field for 
Jc{H,T) measurements at Ec = 10~^ V/m (left graph) and Jc{H,T,e) measurements 
at i?c = 5 X 10"'' V/m (right graph). 



Fp{H,T,e)/Fp^,^{T,e) in: 

FAH,T,e) ^ h'-' (1 - hf ^ h^-' (1 - hf 

i^pmax(T,e) 0.20-5(1-0.2)' 0.286 ' ^ ^ 

where h = H/H*2{T,e). The resulting normahzed bulk pinning force dependencies on 
reduced field are shown in figure |21 

The left graph shows the dependence for the Jc{H,T) measurements at all 
temperatures (and constant strain). The right graph shows the reduced dependencies 
for the Jc{H,T,e) for all temperature and strain values. The fit values for H*2{T) 
and -Fpmax(^) for the Jc{H,T) measurements are summarized in figure El including 
a calculated dependence for H*2{T) using the Maki-De Gennes relation p2j) with 
HoH^^iO) = 30.1 T and T*{0) = 16.7 K. It should be noted that these are at -0.15% 
axial pre-compression, as will be shown in Section 0] The high value for if*2(0) in 
comparison to a measured maximum detectable value /io-f^c2(0) = 29.3 T [60] can be 
explained by inhomogeneities causing a downward tail in Kramer plots thereby resulting 
in an overestimate of the extrapolated H*2{0). It should be emphasized that this high 
value is required to accurately describe the measurements in the range h < 0.8. 

The bulk pinning force curves in figure El demonstrate that p = 0.5 and q = 2 indeed 
yield the correct magnetic field dependence for the entire investigated temperature 
and strain regime. Above h = 0.8 deviations can occur, but these are attributed 
to inhomogeneities. It should be emphasized that the observed correspondence with 
p = 0.5 and g = 2 is not specific for this wire, but holds for all our earlier investigated 
sample materials and also for 'normally' reacted PIT wires as demonstrated by Cooley 
et al. [80]. Note also that the regime above h = 0.8 can be accounted for in the fit 
by using deviating values for p and q. This will obviously lead to a different value for 
the scaling field H*2 and scaling, as demonstrated by a normalized pinning function. 
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Figure 3. Temperature dependence of the effective upper critical field and maximum 
bulk pinning force as measured on the ITER barrel at Ec — 10^^ V/m. 



is still possible as was shown by similar analysis on the same wire by Hampshire and 
co-authors [108]. This will improve the accuracy of the fit above h = 0.8 at the cost of 
generality of the model. However, since this regime involves very low current densities 
that are of negligible importance to the applications, we choose to retain a general model 
at the cost of low current density fit accuracy. 

It is thus concluded that the field dependence is accurately described by p = 0.5 
and g = 2 for all magnetic fields in the range 0.05 < h < 0.8 for all temperatures and 
strains. It is also observed that the Maki-De Gennes relation indeed yields an accurate 
description for the temperature dependence of the effective upper critical field resulting 
from transport critical current measurements. 



3.2. Temperature dependence 

The temperature dependence of the bulk pinning force can now be analyzed. The 
temperature dependencies for the upper critical field and the GL parameter are known. 
This means that only the temperature dependence of the maximum bulk pinning force 
has to be investigated to determine the overall temperature dependence of the bulk 
pinning force. The temperature dependence of the maximum bulk pinning force is 
determined by the values for u and 7 in 



3.2.1. Incorrect temperature dependence of the Kramer /Summers form In the 
Kramer /Summers form z/ = 2.5 and 7 = 2 in (jHUj) and leading to: 

F,[H,T)=C,^^^^h-^[l-h)\ (40) 
kit) 

where Ci = C[^qH*2{0)]^'^ / Ki{Oy . The maximum pinning force is thus given by: 

Fp^ax(T) = 0.286ci^^^i^%P^ 

(T) 
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Temperature [ K ] 



Figure 4. Temperature dependence of the GL parameter and maximum bulk pinning 
force at Ec = 10~^ V/m. The points are derived from the Jc{H,T) characterizations 
and the Hnes are the calculated dependencies using ((TB|l (MDG) and For 
completeness, also the Summers empirical form for the GL parameter is included. 

Ki (0) 

(41) 

There is unfortunately no direct way to extract either C or ki(0) separately from the 
Jc{H,T) results. Comparison to the measurements delivers Ci, H*2{0), T*{0), H*2(T) 
and -Fpmax(^) without free parameters, but ki(0) depends on the value of C. Note 
that C is not a fundamental constant here because Fp is evaluated from the non-copper 
cross-section. 

An overall least squares fit of (pUj) on the measured Jc{H, T) data, using ij,qH*2{0) = 
30.1 T and T*{0) = 16.7 K, resuhs in Ci = 42.8 x 10^ AT/m^ The resulting maximum 
bulk pinning force according to is shown in comparison to the actual data in figurelU 
Significant deviations occur above 9 K, which will also be present in calculated critical 
current densities. 

The deviations in -Fpjnax(^) can be illustrated through the resulting change 
in Ki with temperature. The variation in ki with temperature, according to the 
Kramer /Summers form, can be calculated through (PT|) by using values for -Fpmax(^) 
and H*2(T) from figure El and by inserting a value for ki(0). Choosing C = 
12.8 X 10^ AT^'^m^ as proposed by Kramer (0), results through (jlUj) in ki(0) = 38.6; 
a reasonable value compared to the literature [69]. The resulting values for ni{T), 
calculated from Fp^^xiT) and H*2{T) using pT|l are plotted in figure IH together with 
the microscopic (fT3|) and empirical (fTT|) dependencies. 

The change of i^i{T) of nearly 100% is unrealistic, since a maximum change of 20% 
is expected in the weak, and 50% in the strong coupling limit as discussed in Section 
12.21 This is independent of the choice of C and thus ki(0), since the calculations are 
normalized to ki{0). The combination of = 2.5 and 7 = 2 therefore results in an 
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Figure 5. Kramer plot of measured Jc{H,T) data, including linear fits according 
to and the overall Jc{H,T) fit according to (gOI) at = 10~^ V/m. The graph 
emphasizes the failure of the Kramer/Summers model to account for the correct slope 
at higher temperatures, which leads to significant deviations in the calculated Jc{H, T) 
values. 



unrealistically large change in fi;i(T), calculated from F^^a_^{T) and H*2{T). 

The Kramer function is given by ((Zj) and the slope is therefore C/ki{T). The large 
change in ki with temperature thus translates to a larger change in the slope of a Kramer 
plot at higher temperatures. This is visualized by the Kramer plot given in figure 

The overall Jc{H,T) least squares fit is optimized at 4.2 K in this graph, to 
emphasize the deviations that occur at higher temperatures. The error in the overall 
temperature dependence in the Kramer /Summers model, through the use ofu = 2.5 and 
7 = 2, causes significant deviations in the slope of a Kramer plot at higher temperatures 
and thus in large errors in the calculated Jc{H,T) values. 

Deviations in the Kramer fit were observed also in our earlier analysis [109, 110] 
and were attributed to incorrect (empirical) temperature dependencies of the effective 
upper critical field and the GL parameter. Now that H*2{T) and ki{T) are replaced 
by microscopic based alternatives, it has to be concluded that the observable errors in 
FprnaxiT) stcm from incorrect values for u and 7. This is supported by the unrealistic 
large change in the calculated ki with temperature, surpassing the strong coupling limit. 



3.2.2. Improved temperature dependence In the above analysis it became clear that the 
temperature dependence of the bulk pinning force using u = 2.5 and 7 = 2 is incorrect. 
It is therefore required to determine the optimal values for both constants that yield the 
highest accuracy for the calculation of the temperature dependence of the maximum 
bulk pinning force. The Fp^g^^(T) data are, in the general form, described by: 



Fpmax(T) = 0.286Ci 
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Figure 6. Combinations of v and 7 and the overall least squares fit error in the 
maximum bulk pinning force. 



1/^0^02 (0)] 

«1 (0)^ 



u 



0.286C 



MDG {ty-^ (1 - t^) 



(42) 



Relation (jl^ is fitted against the -Fpmax(r) data from Section 13.11 Values of u are 
chosen around 2.5 and the value for 7 is least squares fitted. Various combinations of 
z/ and 7 are able to describe the measured Fp^s^^iT), each with a specific remaining 
overall least squares error, as plotted in figure El The possible combinations of u and 7 
are related by 7 = 0.880z/ — 0.756. The overall error shows a distinct minimum of nearly 
zero for the combination u = 2 and 7 = 1 and these values therefore yield optimum 
accuracy for the temperature dependence of the maximum bulk pining force. Relation 
(1121) 5 using z/ = 2 and 7 = 1 is plotted together with the -Fpmax(7') data from Section IHTTl 
in figure [71 The accuracy in the calculated -Fpmax(^) is significantly improved compared 
to the use of z/ = 2.5 and 7 = 2 in figure El 

The temperature dependence of the GL parameter, resulting from Fpmax{T) and 
H*2{T) using with z/ = 2 and 7 = 1 can also be calculated. The relation between C 
and /«i(0) now differs from the Kramer form (0). The value for Ki(0) is chosen consistent 
with the previous analysis, i.e. ki (0) = 38.6, leading to C = 1.81 x 10^ AT^'^m^. This 
calculation yields a ki (T) as plotted in figure [71 The GL parameter now changes by 
about 35% which is in agreement with expectations (i.e. between 20 and 50%, see Section 
I2.2.3|) and, moreover, follows the microscopic based description. The overall accuracy 
of the temperature dependence can again also be demonstrated using a Kramer plot. 
Using u = 2 and 7 = 1 in ()36|1 results in an overall description: 



Ki (0) k (t) ^ 

CiMDG (t) (1-t') /i°-^(l-/i)'. 




(43) 
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Figure 7. Improved temperature dependence of the GL parameter and maximum 
bulk pinning force at Ec = 10~^ V/m. The points are derived from Jc{H,T) 
characterizations and the hnes are the calculated dependencies using ()15|) (MDG) and 
(I42|l with v = 2 and 7=1. For completeness, also the Summers empirical form for 
the GL parameter Hll() is included. 



and the linear Kramer function becomes: 

/k (H, T) = ,^^o.,^Z. ,^,10.25 /^0 (H:, (T) - H) . (44) 
Ki (T) [fioH*^ (T)] 

Similar analysis as in the previous Section yields from the overall Jc{H, T) least squares 
fit Ci = 42.6 X 10^ AT/m^. The small difference in Ci compared to the Kramer /Summers 
form (42.8 x 10^ AT/m^) stems from the fact that in figure El the fit was optimized for 
the 4.2 K results, whereas now the fit is optimized over the entire temperature range. 
The improved functionality yields slopes that are consistent with the measured results 
over the entire temperature range as is seen in figure |H1 The only significant deviations 
occur close to H*2{T), the regime where non-linearities in Kramer plots originate from 
A15 inhomogeneities. 



3.3. Strain dependence 

Now that the magnetic field and temperature dependence are established, the strain 
dependence is introduced by stating that the effective upper critical field, the critical 
temperature and the GL parameter depend on strain. This leads through (j35|l to: 

Ki (T, e) 

where H*,{T,e) = i/*2^(T)s(e), r,*(e) = T:^s{tf/^ and «:i(r,e) = «:i^(T)s(e)". It is 
assumed that w = ?> since this was shown earlier to be consistent with measurements, 
and a has to be determined from experiment. In normalized form becomes [through 

Fp {H, T, e) = Cis (e)'"" MDG {t) (l - t^) h^'^ (1 - hf , (46) 
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Figure 8. Kramer plot of measured Jc{H,T) data at Ef. = 10^'"' V/m, including 
linear fits according to (|44|l and the overall Jc {H, T) fit according to (|43(l . The 
improved overall temperature dependence using v = 2 and 7 = 1 results in an accurate 
description across the entire temperature range. 



where Ci = C[/ioi^c2m(0)]V'«im(0). 

First the function for s(e) will be determined. It will initially be assumed that 
a = 1, since this is consistent with Ekin's [71] statement that Fpinax(4.2 K) oc H*2{e), 
as well as our earlier analysis. This choice for a has to be confirmed with experimental 
results. The improved and re-normalized deviatoric strain model is chosen as description 
for s(e), i.e. relation (j22I)- An accurate Jc{H,T,e) measurement with a high strain 
resolution is required to determine the parameters in (j22I), i-e. Cai, Ca2, eo,a and 
em- The mechanical parameters are therefore determined from a measurement on 
the Furukawa wire at 12 T and 4.2 K, which are measured with the required strain 
resolution [107], as shown in figure IHl The points are the Jc(e) results at 12 T and 4.2 
K at a criterion of = 10~^ V/m. The line is calculated from ()46|) combined with (j^^ 
using the mechanical parameters depicted in the graph, « = 1, Ci = 47.8 kAT/mm^, 
/ioi^*2m(0) = 30.7 T and T*^{0) = 16.8 K. 

It should be emphasized that Ci, H*2^{0) and T*^{0) are solely used as fit 
parameters to scale Jc- Since no field or temperature dependent measurements were 
performed on this specific sample they are arbitrary, but not unique, parameters 
resulting from a least squares fit of the 12 T, 4.2 K data points. As mentioned this 
sample was not heat treated together with the U-spring and barrel samples, rendering 
differences of a few percent in Ci, H*2^{0) and T*^{0) possible. In addition, the applied 
criterion of = 10~^ V/m is higher than used for the Jc{H,T) characterizations 
{Ec = 10^^ V/m), translating to slightly higher values for the fit parameters, since 
effectively Sn richer A15 sections will be probed. 

Figure M shows that the improved deviatoric strain model does accurately account 
for the measured behavior. In contrast to the original deviatoric strain model, it is able 
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Figure 9. Axial strain dependence of the critical current density at — 10^^ V/m, 
12 T and 4.2 K measured on a circular bending spring and calculated from (|46|l and 
(|22|l using parameters as depicted in the graph. 



to fit the entire strain curve and not just the compressive part. Note that the term 
Ca2 effectively tilts the calculated curve around H*2^{T) at = 0. This causes the 
fit values for Cai and eo,a to change compared to the earlier model. Overall, however, 
(HD) accurately describes the strain dependence and a = 1 in (jlUj) appears to yield the 
correct strain dependence of the bulk pinning force. 

Now that s(e) is confirmed to be accurate, -Fpmax(e) oc s{ey~'^ has to be investigated 
to find the value for a from experiment, which can be done using earlier U-spring 
results. Note that experimental limitations of the U-spring only allow for a relatively 
high criterion of = 5 x 10~^ V/m to be considered reliably [107]. The power 2 — a 
is determined by plotting the -Fpmax(e) values from Section ITTl versus axial strain and 
comparing to s(e)^~". This is demonstrated in figure ^1 for 4.2 K results, yielding 
a = 1.0 ±0.1, in agreement with Ekin's result. Similar analyses for higher temperatures 
suggest that a might decrease with temperature. At higher temperatures, however, 
the -Fpmax(e) analyses become significantly less accurate since they are based on much 
smaller current densities and limited data. Also inhomogeneity effects, which manifest 
themselves in the investigated magnetic field range at higher temperatures, complicate a 
more precise determination of 2 — a. The results are therefore not sufficiently accurate to 
come to unambiguous conclusions on a possible temperature dependence of a. Moreover, 
strain dependency results are, for now, only available for inhomogeneous samples, 
rendering further speculations not very useful. With the available experimental results 
it can be assumed that a = 1 for all temperatures. 

It is concluded that strain dependence can be introduced in the bulk pinning force 
using a single dependence -Fpmax oc H*2(T, e)^, provided that a = 1 for all temperatures. 
This is in contrast to the usual conclusion that H*2{T) and H*2{e) have different 
dependencies. This single function is enabled through the use of values for u and 7 
that accurately agree with measured results and by including the strain dependence of 
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Figure 10. Maximum bulk pinning force as function of applied axial strain at 
i?c = 5 X 10"'* V/m and 4.2 K, measured on the U-spring. 



the GL parameter. Figure ^1 implies that ki has a similar strain dependence as H*2. 
4. Scaling of measured results 

In the previous Sections an overall scaling function for Fp{H,T,e) was postulated 
and confirmed by systematic analysis of the separate dependencies on magnetic field, 
temperature and strain. The proposed function is given by (j46p . with a = 1. Solving 
the implicit Maki-De Gennes relation can be avoided by using its approximation 
MDG{t) = (1 — t^-^^) which is sufficiently accurate for practical applications. Using 
this approximation changes the overall function to the simple form: 

Jc {H, T, e) = -^s (e) (l - f-^') {l - f) h'-' {I - h)\ 
with 

Ci = C[/Xoi/:2m(0)]7'^lm(0), 

t^T/T;(e), h^H/H:,{T,e), 
H:,{T,e) = H:,^ (O).(e) {l - t'-'') , 

T: (e) = T:^s ie)K (47) 

in which s(e) is given by p2|l . The superconducting parameters Ci, H*2^{0) and T*^{0) 
and the deformation related parameters Cai, 0^2, eo,a and em have to be determined 
experimentally. 

The new scaling relation can be applied to the entire existing Furukawa wire 
Jc{H,T,e) data-set at = 5 x 10^^ V/m which results in the mechanical and 
superconducting parameters given in tabled The axial thermal pre-compression (cm) of 
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Table 1. Parameters for the calculation of Jc{H,T,e) at i?c = 5 x 10 V/ra. 



Deformation 


Superconducting 




C'al C'a2 Co, a 




Ci 


47.6 6.4 0.273 


30.7 16.8 


46.3 



Table 2. Axial thermal pre-compression in percent at T = 4.2 K for a Furukawa wire 
mounted on an experimental setup. 

ITER barrel U-spring Pacman 

-0.153 -0.045 -0.054 



the A15 volumes in the Furukawa wire when mounted on the specific experimental setups 
are given in table |21 These axial pre-compression values are low compared to a wire that 
is not mounted on a substrate (for which the pre-compression usually is between —0.2 
and —0.4%). This is a result of the low thermal contraction of the substrate material 
(Ti-6A1-4V) which determines the overall thermal contraction after wire mounting and 
which matches the thermal contraction of NbaSn. Differences between the thermal pre- 
compression on the ITER barrel and the strain devices are attributed to the differences 
in mounting temperature of the wire (epoxy hardening at room temperature for the 
ITER barrel versus soldering at about 200°C for the strain devices). 

5. Discussion 

5.1. Scaling accuracy 

Through systematic comparison with measured results deficiencies in existing scaling 
relations can be highlighted. Improvements can, however, only be made empirically 
through a change of the bulk pinning force dependency on both H*2{T) and ki(T) 
simultaneously. The resulting relation describes the measurements accurately (also for 
a multitude of other wires) [4,9]. The values z/ = 2 and 7 = 1 are similar to those 
found earlier [4,9,10] through minimization of the error in the overall Jc{H,T,e) least 
squares fit. In the earlier analysis it was assumed that the deviating values for u and 
7 originated from a wrong temperature dependence of Hc2 and ki. It is now concluded 
that the Kramer form, which results from the Labusch [20] function for Cqq, is in 
fact not accurate. These theories have been questioned in the literature, specifically 
by Brandt [21] and Dew-Hughes [17], emphasizing the difficulties in understanding 
and modeling the flux-line to pinning site interactions. A reasonable temperature 
dependence description can be achieved for the combinations 7 = O.SSOu — 0.756, but 
the overall error shows a distinct minimum at z/ = 2 and 7 = 1. For a proper connection 
to theory these values should result from a fundamental form of these interactions. At 
present, to our knowledge, such a model is not available. 
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The validity of the new scahng relation was tested on a database derived from four 
samples, one of which was heat treated separately. The samples were characterized in 
three experimental setups. Deriving a singular parameter set for the Furukawa wire is 
complicated by small differences between the samples and characterization methods. A 
better fit to the U-spring strain results alone can, for example, be found by the use 
of slightly different parameters. This stems from small, for now unclear, differences 
between the separate samples and experiments. Ideally all measurements should be 
performed on a single instrument and on a single wire and at one single, low criterion. 
The latter is preferred to remove inaccuracies arising from parasitic current subtraction 
and for consistency with the more standard criterion of Ec = 10~^ V/m. Nevertheless, 
the accuracy of a very similar function (using the earlier deviatoric strain model) was 
investigated extensively in previous work and it was concluded that a standard deviation 
below 2% can be obtained [2]. The applicability of the proposed scaling function was 
validated for Internal-Tin wires during the ITER benchmark tests in our earlier work [9], 
and in unpublished analysis of Powder-in- Tube wires. 

5.2. Minimum required dataset 

With a scaling relation available the question arises as to what is the minimum data-set 
that is needed to retrieve the entire critical surface. For the temperature dependence 
two H*2(T) points are required to determine the entire H*2(T) dependence through the 
MDG relation. These should be determined by transport Jc measurements since these 
yield the correct H*2 and T* required for scaling. From figure |H1 it is seen that preferably 
at both ends of the temperature range the slope of the Kramer plot has to be determined 
to extrapolate to H*2{T). This is most simply done by a Jc{H) measurement at 4.2 K at 
a few magnetic field values to check for linearity and thus validate the choice of g = 2. 
In addition this should be done at a higher temperature, e.g. 12 K might be a good 
compromise between temperature range and available magnetic field in a laboratory 
magnet. Non-linearities above h = 0.8 should be neglected since they most probably 
result from inhomogeneities. The results can then be extrapolated linearly in a Kramer 
plot to yield H*2{4:.2 K) and H*2{12 K). The MDG relation (or its approximation) then 
yields H*2{0,e = em) and T*{0,e = em)- Jc{H,T,e = e^) and Ci are known through 
(jTTj) using a least squares fit to the Jc data-set. 

The strain dependence can then be determined by a 4.2 K strain measurement with 
sufficient strain resolution to accurately determine the mechanical parameters (e.g. as 
in figure IH)) . The magnetic field for the strain characterization should preferably be 
chosen such, that over the entire strain range of interest, the critical current remains 
above about 10 A to avoid inaccuracies at low current measurements. A comparison to 
the Jc(iJ, 4.2 K) results above will yield em and H*2^{0) and T*^{0) can be calculated. 
The entire critical surface up to the magnetic field where inhomogeneity effects start 
to appear (i.e. h = 0.8) is then known. The characterizations should preferably be 
done using one sample and one instrument to avoid experimental errors. The required 
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data-set can thus be summarized as follows: 

(i) An Ic{H) measurement at 4.2 K from about 5 T up to the maximum available 
magnetic field (see figure IH}, preferably 12 to 15 T. 

(ii) An Ic{H) measurement at about 12 K from 1 T or as low as possible, up to the 
maximum available magnetic field. The selectable temperature depends on the 
available magnetic field range (see figure |H1). 

(iii) An /c(e) measurement at 4.2 K and at a magnetic field where the critical current 
remains significant (e.g. > 10 A) over the entire strain range of interest (see figure 

6. Conclusions 

A new, generally valid scaling relation is proposed which is based on a Kramer type 
pinning description. A number of improvements are made in comparison to the Summers 
form. A first improvement is a replacement of the empirical temperature dependency 
for the upper critical field with a microscopic based alternative. A second improvement 
is to replace the empirical temperature dependence of the Ginzburg-Landau parameter 
by a microscopic based alternative. A third improvement is made by replacing the 
power-law strain dependency with a new empirical deviatoric strain description that 
accounts for the effects of the third strain invariant and thus for asymmetry in the 
axial stain dependency curve. The use of an explicit normalized strain dependence 
term s(e) enables straightforward implementation of optional future, more fundamental 
based, strain descriptions in the general scaling relation. A fourth improvement involves 
accounting for the inhomogeneities that are present in the A15 regions in wires by 
recognizing that the critical current scales with a bulk, inhomogeneity averaged critical 
field, representing some average of the upper critical field distribution that is always 
present in the A15 in wires. 

A systematic comparison of measured critical current as function of magnetic field, 
temperature and strain with scaling relations shows that fiux-line shear can be assumed 
to be the main de-pinning mechanism and that tails in Kramer plots can reasonably 
be attributed to A15 inhomogeneities. This results in a generally valid magnetic field 
dependence of the Kramer form, i.e. Fp oc h^'^{l — hY within the range 0.05 < h < 0.8. 

Through the analysis of the magnetic field and temperature dependence of the bulk 
pinning force it is shown that the Kramer/ Summers form, i.e. Fp oc [if*2(T)]^'^/Ki(T)^, 
fails to describe the observed temperature dependence with sufficient accuracy. The 
earlier suggestion that this followed from an empirical temperature dependence of 
the upper critical field and the GL constant was contradicted through the use of 
microscopic based relations for both. The change in ki{T) in the Kramer /Summers 
form is significantly larger than 50% as is expected from microscopic theory in the 
strong coupling limit. The improved temperature dependence of the form Fp oc 
[if*2(T)]^//ti(T) does accurately describe the measurements. The resulting change in 
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K'i{T) follows the microscopic form and is within the strong coupling limit. The new 
dependency is in disagreement with existing relations for the flux-line shear modulus 
6*66 • It is suggested that the temperature dependence of Cee as given by Labusch needs 
to be reevaluated. 

The new temperature dependence of the bulk pinning force allows for an 
unambiguous introduction of the strain dependence of the bulk pinning force, without 
the need to include separate functions to describe the strain dependent upper critical 
field and temperature dependent upper critical field, as is usually done in the literature. 
The measurements indicate that the Ginzburg-Landau parameter has similar strain 
dependence as the upper critical field. 

The improved deviatoric strain model, which includes asymmetry, yields an 
accurate description of measured behavior. Though empirical, it is based on conclusions 
resulting from more fundamental approaches and it recognizes the three dimensional 
nature of strain. 
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